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ABSTRACT 

A nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants is employed to 
re-examine the nonthermal properties of the remnant of SN 1987A for an extended evolutionary period 
of 5-100 yr. It is shown that an efficient production of nuclear CRs leads to a strong modification of 
the outer supernova remnant shock and to a large downstream magnetic field « 20 mG. The shock 
modification and the strong field are required to yield the steep radio emission spectrum observed, 
as well as to considerable synchrotron cooling of high energy electrons which diminishes their X- 
ray synchrotron flux. These features are also consistent with the existing X-ray observations. The 



expected 7-ray energy flux at TeV-energies at the current epoch is nearly e 1 F 1 



4 x 10 



-13 



erg 



cm s under reasonable assumptions about the overall magnetic field topology and the turbulent 
perturbations of this field. The general nonthermal strength of the source is expected to increase 
roughly by a factor of two over the next 15 to 20 yrs; thereafter it should decrease with time in a 
secular form. 

Subject headings: acceleration of particles — ISM: individual (SN 1987A) — supernova remnants — 
X-rays: individual (SN 1987A) — gamma rays: observations 



1. INTRODUCTION 

The Supernova (SN) which occurred in 1987 in the 
nearby Large Magellanic Cloud was the first object of its 
kind in modern times whose evolution could be spatially 
resolved as a function of time; this includes the char- 
acterization of the progenitor star. SN 1987A has been 
extensively studied in all wavelengths from the radio to 
the 7-ray ran ge. For a summ ary of its global character- 
istics see e.g. iMcCravl (|1993l) . 

The properties and the future evolution of this super- 
nova can be viewed from different physics points of view. 
The present work concentrates on the nonthermal char- 
acteristics and, in particular, on the particle acceleration 
aspects and the e xpected 7-ray emission. It extend s two 
ea rlier studies by iBerezhko Sz Ksenofontovl ((2000), and 
bv lBerezhko fc Ksenofontov! (|2006h . hereafter referred to 
as BK, to include the subsequent Ch andra observation s 
of the soft and hard X-ray emission (jParket al.ll200l . 
as well as the most recent and very de tailed observa- 
tions of the radio continuum emission bv IZanardo et all 
(2010). At the same time a qualitative and in fact semi- 
quantitative prediction of the future nonthermal emission 
from this source for the next decades is attempted. 

The radio synchrotron emission and its spectral index 
are important for the determination of the nonthermal 
energy density in the source. In particular the softness 
of the observed radio synchrotron spectrum shows that 
the radio electrons are accelerated in a rather weak shock. 
This is the plasma subshock formed when the accelerated 
nuclei modify the supernova remnant (SNR) blast wave 
by their pressure gradient, so that the main gas compres- 
sion occurs in an extended CR precursor that does not 
accelerate radio electrons. This nonlinear modification 
determines the injection rate of the nuclear particles and 
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vize versa (see BK and section 3). 

The new Chandra data in addition suggest that the 
hard X-ray emission is predominantly nonthermal. Since 
at the corresponding electron energies synchrotron losses 
are of overriding importance for the shaping of the X-ray 
synchrotron spectrum these new observations also deter- 
min e the effective m agnetic field strength in the source 
(see Berezhko 2008, for a review). It turns out that these 
recent observations can be consistently explained with 
only minor modifications of the previous theoretical con- 
siderations. In particular in can be shown that the cir- 
cumstellar environment will lead to a further increase of 
the expected 7-ray emission in the next decades by a 
modest factor of only about 2, before the source will go 
into a steady emission decrease as the SNR shock leaves 
the immediate circumstellar environment and propagates 
secularly into the progenitor star's unperturbed Red Su- 
pergiant (RSG) wind region. 

In comparison with the other well-studied young, re- 
solved SNR in the early sweep-up phase, which also ex- 
hibits nonthermal X-ray emission and a steep radio syn- 
chrotron spectrum, the Galactic object G1.9 -I- 0.3 that 
is presumably the remnant of a type la SN in an es- 
sentially uniform interstellar medium (ISM), SN 1987A 
has a complex immediate circumstellar environment due 
to the interactions of its latest wind phases and the 
stellar radiation field before the explosion. This leads 
to a very non-uniform density structure interior to the 
stellar wind bubble formed by the progenitor in its 
RSG phase. Therefore G1.9 + 0.3 and SN 1987A pre- 
sumably constitute the extremes of spatially resolved 
pre-supernova circumstellar environments when studying 
the particle acceleration properties of extremely young 
SNRs. To this extent the present work is complemen- 
tary t o the analysis of the nont hermal properties of G1.9 
+ 0.3 (Ks enofontov et aDl2010f ). As one might expect on 
physics grounds, it turns out that in both cases the SN 
blast wave is already strongly modified by the accelerated 
particles - mainly atomic nuclei, called here cosmic rays 
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(CRs) — from the very beginning of SNR formation, al- 
though the total amount of nonthermal energy, like that 
of thermal energy, is still only a percent fraction of the 
total explosion energy. This shock modification is asso- 
ciated with a strong amplification of the magnetic field 
strength in the swept-up circumstellar material. These 
conclusions are derived from the following observations: 
(i) the existence of nonthermal X-ray emission (see sec- 
tion 4), (ii) the rapi d increase of the radio emission (see 
iMurphv et al.l 120081 for a recent analysis of the radio 
light curve for G1.9 + 0.3), a nd (iii) the anomalously so ft 
radio emission spectra fBK. iKsenofontov et al.1 ([2010)). 
In both cases these evidences come from the observed 
properties of synchrotron emission produced by the elec- 
tron CR component. The electron component therefore 
plays a fundamental role, permitting indirect conclusions 
about key properties of the nuclear energetic particle 
component, like its pressure and its amplification of the 
magnetic field. Given the diffuse radiation field, the 
inverse Compton gamma-ray emission by the electrons 
is then also determined, and a measurement of the to- 
tal gamma-ray emission yields the hadronic gamma-ray 
component. 

SN 1987A and Gl. 9+0.3 are nevertheless distinct 
from other young historic Galactic SNRs like the core 
collapse object RX J1713.7 + 3946 and the type la 
object SN 1006, respectively, in that the latter SNRs 
have already swept up an amount of circumstellar 
matter that is approximately equal to the ejected 
mass. Therefore they are expected by a number of 
authors to already contain an amount of nonthermal 
particle energy that corresponds to s ome 10 % of the 
total h ydrodynamic explosion energy jA haronian et all 
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As a consequence of the softer effective equation of 
state of the combined gas-CR system the dynamical evo- 
lution in this phase is different from that of a purely gas 
dynamic system. Also the increasing dissipation of hy- 
dromagnetic fluctuation fields will presumably influence 
the system (jPtuskin k Zirakashvila l2003i I2005D . not to 
speak of an early onset of radiative gas cooling as a re- 
sult of the increased gas compression in the interio r that 
is hig her than the purely gas dynamic ratio of 4 (Dorfi 
[1991 . 

The observed shape of the remnant, especially at ra- 
dio frequencies is not very far from spherically symmet- 
ric, even though the o bserved emission v aries around the 
remnant's periphery (|Potter et al.l [200 9). This suggests 
that a spherically symmetric model is a reasonable first 
order approximation for a study of the general nonther- 
mal properties of the remnant. 

In the next section an approximate spherically sym- 
metric model of the pre-explosion environment of 
SN 1987A is formulated. This is combined in section 
3 with a theoretical, spherically symmetric, nonlinear 
model for the nonthermal evolution, using the earlier 
work of BK, and with the recent observations of the 
synchrotron emission, mentioned before. The results are 
presented in section 4 and summarized in section 5. 

2. APPROXIMATE SPHERICALLY SYMMETRIC MODEL 
OF THE CIRCUMSTELLAR ENVIRONMENT SNR 1987A 



The initial sh ort outburst of radio emission 
([Turtle et al.l [T987) has been attributed to the syn- 
chrotron emission of electrons accelerated by the SN 
shock propagating in the free wind of the blue su- 
pergi ant (BSG) progenitor star ([Chevalier k Franssonl 
Il987f l . After about 3 years radio emission was detecte d 
again (|Stavelev-Smith et al.l [19921 iGaensler et al.lll997l) . 
as well as a monotonically increasing X-ray emission 
([Gorenstein et al.1 119941 : iHasinger et ail 11996( 1 . This 
second increase of emission has then been attributed 
to the entrance of the outer SN shock into the wind 
bubble of the BSG of density pe, thermalized in a ter- 
mina tion shock at a radial distance Rt =3.1 x 10 17 cm 
(e.g. IBerezhko k Ksenofontovl [20001 . and subsequently 
into an H II region of much denser matter, consist- 
ing of the swept-up wind of a RSG precursor phase 
(|Chevalier k Dwarkadasl [l995(l . The density pr of the 
swept-up RSG wind region is large compared to pe- The 
contact discontinuity between the two winds is assumed 
at Rc = 5 x 10 17 cm, the scale lc of a modeled smooth 
transition is assumed as lc = 0.05i?c- 

These parameters are consistent with canonical values 
of stellar ejecta mass M Ci = 10M©, distance d — 50 kpc, 
and hydrodynamic explosion energy E sn = 1.5 x 10 51 erg 
(e.g. lMcCrav| [l993). During an initial period the ejecta 
material has a broad distribution in velocity v. The 
fastest part of this distribution can be described by a 
power law dM e -Jdv oc v 2 ~ k . The present modeling uses a 
value k = 8.6 appropriate for SN 1987A ([McCravfl 19931) . 
The general picture is that the interaction of the ejecta 
with the circumstellar medium (CSM) creates a strong 
shock there which heats the thermal gas and accelerates 
particles diffusively to a nonthermal CR component of 
comparable en ergy density. 

According to I Chevalier k Dwarkadasl (|1995l) the CSM 
at r > Rc includes three subsequent regions: the H II 
region, an equatorial ring and, beyond the ring, a free 
RSG wind (followed further out by a rarefied main se- 
quence wind bubble which is of no concern here). Nei- 
ther the contact discontinuity between the two winds is in 
reality a spherical surface at radius Rc, nor is the equa- 
torial ring a sph erical shell, as the optical observations 
clear ly indicate ([McCravl 119931 : IChevalier k Dwarkadasl 
1995). Nevertheless, the present nonthermal model ap- 
proximates the mass distribution per within the ring by 
a spherically symmetric shell 

Per = p m exp[-(r - # R ) 2 /Z R ], (1) 

where p m « M s h/(47T 3 / 2 i? 2 h ZR) is the central (maximal) 
density of the shell; M s h, -Rr and Zr represent the total 
mass, the radius and the width of the shell, respectively. 
Below, the values -Rr = 7 x 10 17 cm and Zr = 0.25-Rr 
are used. Approximately Zr corresponds to the interval 
of radial distances where the dense cool shell, associated 
with the equatoria l ring, is positioned around the source 
(|Crots et al.l[l995l ). 

In the same spirit, the density profile of the RSG mat- 
ter consists of three corresponding terms: 

PR = p\\H(R K - r) + per + p w H(r - R R ). (2) 
Here H(x) is the Heaviside function. 

This results in the following model for the CSM density 
distribution at the distances r > i?x in the form 

t r Pb + PR Pb - pn , , r - R c , . 
pVo = g 2 ~l ' 
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For the gas number density N g = p/rn p the following 
values are adopted: = 0.29 cm' 3 , = 280 cm~ 3 , 
where m p is the mass of proton. These values, as was 
shown by BK, provide a good compromise between the 
SN shock dynamics seen in the radio and X-ray emis- 
sions. 

Since the density of the free RSG wind can be ex- 
pressed in the form 



M 



(4) 



where M is the RSG mass loss rate and V w denotes the 
RSG wind speed, the mass of the shell can be expressed 
in the form 



M„i 



MRr 47rpii(g| ~ R 3 C ) 



(5) 



It is assumed that the center of the spherically smoothed 
ring is positioned exactly at the boundary r — i?R be- 
tween the HII r egion and the free RSG wind region. 

According to iBlondin k Lundqvistl (|1983h the RSG 
mass loss rate M is highly nonuniform: it is M = 
8 x 10~ 5 Af Q yr" 1 in the equatorial region within ±9°, 
and M = 4 x 1O _6 M0 yr _1 in all other directions. In 
the present spherically symmetric approximation its av- 
erage value is then M = 1.6 x 10~ 5 M Q yr _1 . 

For the conventional value V w — 5 km s _1 of the RSG 
wind speed the masses of the shell is M s h =0.5x M . 

3. PARTICLE ACCELERATION MODEL 

To describe the observed properties of nonthermal 
emission of SN 1987A a nonlinear kinetic theory is used 
here. It couples the particle acceler ation process with the 
hydrodynamics of the t hermal gas (|Berezhko et alj[l996t 
iBerezhko k V olk 2000). Therefore, in a spherically sym- 
metric approach it is able to predict the evolution of gas 
density, pressure, mass velocity, together with the en- 
ergy spectrum and the spatial distribution of CR nuclei 
and electrons at any given evolutionary epoch t, and the 
properties of the nonthermal radiation produced in SNRs 
due to these accelerated particles. 

The model does not contain the reverse shock propa- 
gating in the SN ejecta. The reverse shock is not ex- 
pected to be an efficient CR accelerator, because the 
magnetic field value is expected to be very weak as 
a result of the large expansion of the exploded stel- 
lar material that leads to inefficient high energy CR 
acceleration. Even if one assume equally efficient CR 
injection/acceleration on the forward and the reverse 
shock, the expected energy content of CRs accelerated 
by the forward shock is considerably larger compared 
with the reverse shock except during the transition phase 
from the free expansion to th e Sedov SNR phase (see 
IBerezhko k Ksenofontovl 12000. for more details). Since 
SN 1987A is still far from the Sedov phase one can neglect 
the contribution of the reverse shock to the nonthermal 
emission of the remnant. In addition the analysis of the 
obser ved fine structure of the X-ray and the radio emis- 
sion (|Zhekov et al.l 120101) gives evidence that they are 
produced by the forward shock (blast wave). 

In the particle acceleration model the scattering prop- 
erties of the CSM are required. Given the enormous 
shock velocities in excess of w 5000 km s' 1 (see next 



section), Bohm diffusion appears a most reasonable ap- 
proximation. It was also used by BK. In this limiting 
case the scattering mean free path equals the particle 
gyro radius and is therefore inversely proportional to the 
mean magnetic field strength B. 

A rather high downstream magnetic field strength 
£?d ~ 10 mG is required to reproduce the observed ra- 
dio and X-ray spectra (BK). It is far from obvious that 
the dynamical interaction of the RSG and BSG wind 
systems, responsible for the inferred circumstellar struc- 
ture, can also lead to such high magnetic field strengths. 
It is much more likely that the required strength of the 
magnetic field has to be attributed to nonlinear field am- 
plification at the SN shock by the CR acceleration pro- 
cess itself. Ac cording to plasma p hysic al consider ations 
going back to iLucek k BelH (|2000f ) and US] (|2004f >. the 
existing CSM magnetic field can indeed be significantly 
amplified at a strong shock by CR streaming instabili- 
ties. In fact, for all the thoroughly studied young SNRs, 
the ratio of magnetic field energy density Bq/8it in the 
upstream region of the sh ock precursor to t he CR pres- 
sure P c is about the same (fVolk et al.ll2005f ). Here Bq is 
the far upstream field in the precursor, presumably am- 
plified by the energetically dominant CRs of the highest 
energy. Within an error of about 50 percent the empiri- 
cal relation 



B = y/2n x io- 2 p y s 2 



(6) 



holds, where V s denotes the SN shock speed. This is 
quantitatively close to the above-mentioned field value. 
This is the same field value as used by BK. The factor 
10~ 3 instead of 10~ 2 in Eq.(4) of BK is simply a misprint: 
this is clearly seen from a comparison of Fig.l of BK with 
Fig.l of the present paper. 

For the following it is assumed that this relation also 
holds for an individual young source like SN 1987A dur- 
ing its temporal evolution, where V s (t) depends on the 
source age t. For time scales of the order of the propaga- 
tion time of the shock through th e shock precursor, i.e . 
about 10% of the system age (e.g. IBerezhko et al.|[l996|) 
this appears as a reas onable approximation . 

As reviewed earlier (|VoHd 120041: lBerezhkoll2005l 12 008). 
and elaborated most recently in detail ( Berezhko et al.l 
2009; Ksenofo ntov et all 12010ft . the key parameters of 
the theoretical model (proton injection rate ij and elec- 
tron to proton ratio K cp below the synchrotron cool- 
ing range) can be determined in a semi-empirical way 
from a fit of the nonlinear theoretical solution to the 
observed synchrotron emission spectrum. For the sake 
of simplicity the values of these parameters are usually 
assumed to be constant during SNR evolution. Also 
BK proceeded in this way. The present paper assumes 
such a simplification to be the main reason why a cor- 
responding theoretical model fits the data only on aver- 
age. This concerns the shape of the radio spectrum and 
the time-dependence of th e radio flux over an extended 
period of SNR evolution (|Zanardo et al.l [20101) . In the 
case of SN 1987A, on the other hand, it is natural to ex- 
pect substantial time- variations of these injection param- 
eters since the SN shock propagates through a strongly 
nonuniform CSM whose physical parameters at the shock 
front are expected to change in time. Therefore, in con- 
trast to BK, time dependent injection parameters r](t) 
and K cp (t) will be admitted here. Their values are de- 
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Fig. 1. — Shock radius R B (solid), shock speed V a (solid), gas 
density N e (dash-dotted) and upstream magnetic field Bo (dashed) 
at the current shock position as a function of time since SN explo- 
sion for = 7 X 10 17 cm. The dotted vertical line marks the 
current epoch. The observed radius a nd speed of the SN shock, 
as determined by radio observations (Ne ct al. 2008) are shown 
as well. The scaling values are Ri = i?x = 3.1 X 10 17 cm and 
Vi = 28000 km s" 1 

termined from the fit to the measured synchrotron data. 
The main physical factor which determines the injection 
efficiency r](t) is the structure of the magnetic field up- 
stream of the shock. Injection is expected to be progres- 
sively less efficient when the magnetic field component 
tangential to the shock surface becomes more relevant 
(jVolk et al.ll2003l ). Such a situation is expected in the 
most compressed CSM region, that is, in the ring region. 
It is important to note that a reliable semi-empirical es- 
timate for the values of these parameters for any given 
SNR evolutionary phase is possible if the measured syn- 
chrotron spectra in the radio and X-ray bands are avail- 
able for this particular p hase (e.g. iBerezhko et al J 120091 : 
iKsenofontov et all 12010 ) . For SN 1987A good quality 
measurements of the radio spectrum, as well as of the 
X-ray fluxes in the soft (0.5-2 keV) and the hard (3-10 
keV) energ y ranges exist now fo r the full evolutionary 
period (e.g. IZanardo et al|[201O ). The remaining prob- 
lem is that the nature of the observed X-ray emission of 
SN 1987A is not known unequivocally. In the extreme 
case, when only an upper limit for the nonthermal (syn- 
chrotron) emission is known, the data yield only a lower 
limit for the magnetic field strength Bo(t) and a cor- 
responding value of the electron to proton ratio (BK). 
It is, however, noted here that the time dependence of 
the hard X-ray flux differs from that of the soft X-ray 
flux, and that it is very close to the time dependence of 
the radio emission flux. This fact can be interpreted as 
evidence that the hard X-ray emission is predomi nantly 
of a nonthermal nature (e.g. IZanardo et ; 1 1 .1 [Uo 1 () :■ . The 
considerations below are based on this interpretation. 

The present model calculations start at the SNR evo- 
lutionary epoch t = 1000 d, when the outer SN shock has 
reached a radius Ri = Rt and a speed V\ = 28000 km 
s _1 . These values of R s and V s correspond to the end of 
the SN shock propagation in the free BSG wind region 
r < Rt (BK). The contribution of CRs accelerated in 
the region r < Rt is neglected, because the number of 
CRs produced in the region r > Rt becomes dominant 
very soon, because of the high gas density there. 



4. RESULTS AND DISCUSSION 

The calculated shock radius i? s and shock speed V s , 
shown in FigQ] as a function of time, are in satisfac- 
tory agreement with the values obtained on the basis 
of radio measurements. Up to the year 2030 the shock 
speed decreases due to the increase of the CSM density 
and then, after crossing the maximal CSM density, it in- 
creases again from V s = 3300 km s _1 to V s ~ 4400 km 
s _1 . The reason for this unexpected behavior is that 
during the whole time period under consideration the 
swept-up mass is much lower than the ejecta mass M e j. 
Therefore the main fraction of the explosion energy is 
contained in the freely expanding ejecta. On a small 
time scale the shock is driven by the downstream over- 
pressure: each local relatively sharp increase of the CSM 
density leads to a temporal decrease of the shock speed 
that keeps PqV^ constant. However, on a larger time 
scale the shock is piston driven. 

Note that the actual CSM is essentially non-spherically 
symmetric during the time period considered. Despite of 
this situation the piston-driven shock is expected to be 
roughly spherical. This supports the spherically sym- 
metric approach used here. 

The fit of the theoretical solution to the synchrotron 
spectra, as measu red in the radio and X-ray ranges up 
to the year 2008 (|Zanardo et al.ll20TO) . yields estimates 
for the proton injection rate r)(t) and for the electron- 
to-proton ratio K ep (t). This is shown in Fig|2] The 
fit procedure was de s cribed in detail for sim ilar cases 
(jBerezhko et alll2009l: IKsenofontov et al.ll2010D . There- 
quired proton injection rate r](t) ~ 10~ 3 leads to a signif- 
icant nonlinear modification of the shock: as can be seen 
in Fig(5] the total shock compression ratio a = 5.4 — 6 
is essentially larger, and the subshock compression ra- 
tio a s w 2.8 — 3.6 is lower than the classical value of 4 
for a pure gas shock. Fig|2] also shows that the proton 
injection rate ry changes during the SNR evolution, so 
that it has a local minimum at an age of about 20 yr. 
This is required in order to fit the observed radio emis- 
sion spectra, which are becoming somewhat harder than 
a spectrum calculated with an unchanged injection rate 
- at least up to the present time. During the period from 
t w 13 yr to t rj 30 yr, i.e. beyond the present epoch 
of radio observations, the shock is within the dense shell 
that corresponds to the observed equatorial ring. In this 
spatial region the magnetic field is expected to have an 
enhanced tangential component which should depress the 
nuclear injection rate. 

Thereafter the injection rate is a priori unknown. It 
may increase again to the same level r\ (t) « 10~ 3 . This 
is one of the possibilities calculated below (so-called high 
proton injection rate). However, since the magnetic field 
vector in the free RSG wind is dominated by a compo- 
nent tangential to the SNR shock also the possibility of 
a continuing depressed injection (so-called low proton in- 
jection rate) at a level similar to the value rj(t) k4x 10 -4 
within the ring region is considered. 

The higher injection rate yields for t > 23 yr a corre- 
spondingly larger shock modification, characterized by a 
larger shock compression ratio a and a lower subshock 
compression ratio a s (Fig|5]). 

The strongly modified SNR shock generates a CR spec- 
trum iV oc p -7 , which is very soft at momenta p < m p c, 
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Fig. 2. — Shock compression ratio a (solid), subshock compres- 
sion ratio <j b (solid), proton injection rate rj (dashed), and electron- 
to-proton ratio K cp (dash-dotted) as a function of time. The dotted 
vertical line again marks the current epoch. The thick and the thin 
lines here and in the subsequent figures correspond to an increasing 
(= high) and to a constant (= low) proton injection rate after an 
age of ~ 23 yrs, respectively. 
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Fig. 3. — Synchrotron spectral energy density of SN 1987A, cal- 
culate d for the five evolutio nary epochs (solid curves). The ATCA 
radio {Zanardo ct al. 2010) and dat a for three epochs are shown as 
well, together with Chandra X-ray fPark ct al. 2007) data for two 
epochs (both, in the soft energy range at v = 10 17 Hz and, con- 
nected by dotted lines, in the hard energy range at v = 6 X 10 17 Hz). 
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Fig. 4. — Spectral index of the radio synchrotron emission of 
SN 1987A as a func tion of time, together with the ATCA data 
HZanardo et al.|[20Tot) . 

with index 7 « (<7 B +2)/(cr s — 1) m 2.7. CR electrons with 
such a spectrum produce a radio synchrotron emission 



spectrum S„ oc v~ a with spectral index a = (7 — l)/2 w 
0.9, that corresponds very well to the observations, as 
can be seen in Fig(3J where the synchrotron spectral en- 
ergy density vS v is calculated for five successive epochs 
together with the experimental data. 

The assumed time variation of the proton injection rate 
7/ leads to a time dependence of the radio spectral index 
a as shown in Fig|U and it prov ides a better fit of the ob- 
servations (Zanard o et al.ll201Q[ ) compared with the case 
of constant r\ (BK). In fact the adopted injection rate 
r)(t) at t < 20 yr is the result of the fit to the radio 
synchrotron spectrum. 

For the future, t > 23 yr, the slope of the radio syn- 
chrotron spectrum is expected to be nearly the same if 
the proton injection remains constant at r\ 4 x 10~ 4 
(FigH]). It is expected to become progressively steeper 
due to the decrease of the subshock compression ratio 
a s (t) in case the proton injection rate increases again to 
the level rj « xl0~ 3 . 

The strong downstream magnetic field ~ 20 mG 
leads to synchrotron cooling of the electrons with mo- 
menta p > 10m p c. This makes the high energy part 
of the synchrotron spectrum ( v > 10 13 Hz) very soft 
(see FiglHJ}. At higher frequencies ( v = 10 16 - 10 19 Hz) 
the synchrotron spectrum becomes harder, possibly due 
to a pile- up effect. It hardens the spectrum of acceler- 
ated electrons, that undergo strong synchrotron losses , 
just near its exponential cutoff (e.g. iDrurv et al.lll999h . 
Under this condition the calculated synchrotron flux at 
frequency v w 10 17 Hz, which corresponds to a photon 
energy e 7 = 0.5 keV, is below the measured flux (see 
FiglSJ). This is a required condition because at energies 
e 7 = 0.5-2 keV the X-ray emission of SN 1987A is dom- 
inated by lines and is therefore mainly of thermal origin. 
At higher energies e 7 > 3 keV X-rays are presumably 
of a predominant nonthermal origin. Therefore the fit 
of the measured X-ray flux for e 7 = 3 keV was used in 
the determination of the values of the injection parame- 
ters (r](t) and K cp (t)) and of the amplified magnetic field 
Bo(i) at the beginning of the shock precursor. 

The less modified shock corresponding to a low injec- 
tion rate during the future epoch t > 23 yr yields lower 
cooling of electrons due to the lower downstream mag- 
netic field on account of the lower overall compression 
ratio. This leads to a considerably flatter high-energy 
part of the synchrotron spectrum ( v > 10 13 Hz). 

The calculated 7-ray integral spectral energy flux den- 
sity (SED), shown in FigEl is dominated by the 7r°-decay 
component at all energies. 

It is important to note here that the hadronic SED 
has nevertheless been renormalized by a factor / rc = 0.2 
compared to the amplitude resulting from the spherically 
symmetric model on which the entire calculation is predi- 
cated (e.g. |MOE]iXll2QQ3)- The reason for this a poste- 
riori correction is the fact that the circumstellar environ- 
ment is characterized by a Archimedean spiral topology, 
both in the BSG as well as in the RSG wind bubble en- 
vironment. In particular the radiative cooling and the 
ionization by the stellar radiation field of the progeni- 
tor star will make the entire interaction region of these 
successive wind phases highly turbulent. Therefore there 
will always be regions where the SNR is quasi-parallel, 
that is to say, regions where the shock normal will make 
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Fig. 5. — Integral 7-ray energy flux from S N 1987A, calculated 
for four epochs. The CANGAROO ( arrows) lEnomoto et al.ll2007ft 
and the latest H.E.S.S. (triangles) l IKomin et al .120101) upper limits 
are shown as well. 

an relatively small angle with the local upstream mag- 
netic field direction. Only in these quasi-parallel shock 
regions nuclear particles can be effectively injected into 
the diffusive acceleration process. However, as a conse- 
quence of the overall Archimedean spiral topology, the 
total solid angle of these quasi-parallel shock segments 
will be 47r/ re , where / re is considera bly smaller than 
unity. Following pre vious estimates (|V61k et al.l 120031 : 
iBerezhko et al.l I2003D / re = 0.2 is used in this paper, 
but the actual value of / rc may be even smaller. The 
corresponding uncertainty in the amplitude of the 7-ray 
energy flux remains. 

Since the SN shock is strongly modified the 7-ray 
spectrum at energies e 7 > 0.1 TeV is very hard: F 7 oc 
e 7 ' 8 . At the current epoch the expected 7-ray energy 
flux at TeV-energies is about e 1 F 1 «4x 10 -13 erg cm' 2 
s _1 . and during the next 20 years it is expected to grow 
by a factor of about two with a subsequent temporal de- 
crease F 7 oc R~ x due to the rapid spatial decrease of 
the C SM density p w {R s ) oc R~ 2 (e.g. IBerezhko fc Volkl 
I2000T) . 

At present there exist only upper limits for the TeV 
emission. They have b een obtained b y the CANGAROO 
(lEnomoto et all 12007ft and H.E.S.S. komin et al][20Tol) 
instruments (see FigJ5|). The latest H.E.S.S. upper limit 
is quite close to the SED estimated above. Even with 
the reservations formulated above, this is remarkable and 
justifies in the view of the present authors a deep obser- 
vation of this object. 

The expected time profile F 1 (t) is sensitive to the ra- 
dial profile of the actual CSM density distribution. If the 
dense shell, which in the spherically symmetric model 
represents the matter contained in the equatorial ring, is 
situated at a larger distance, say at i?R = 9 x 10 17 , the 
peak of 7-ray flux is expected to occur 10 years later, 
with an amplitude that is by a factor of 1.3 lower than 
in the former case. 

It is clear from the above consideration that the pro- 
ton injection rate 77 (whose value influences the efficiency 
of CR production) can be estimated from the observed 
shape of the radio emission. Since it is not possible to 
predict the values of 77 for the future evolutionary epochs, 
the presented prediction for the corresponding gamma- 
ray emission is uncertain. Note however that this kind of 



uncertainty is not very big due to the following reasons. 
First of all, it is hard to expect that the shape of the ra- 
dio spectrum will be much more rapidly changing during 
the next 20-30 years than during the previous epochs. 
Therefore the predictions based on the recently deter- 
mined value of 77 should be roughly valid also for this 
period of time. Secondly, the shape of the CR spectrum 
is strongly sensitive to the injection rate at low energies 
e < m p c 2 . Since the radio synchrotron emission is pro- 
duced by electrons with such energies it is sensitive to 
the expected injection rate (Fig |3l4[) . The high-energy 
part of the CR spectrum (e ^ m p c 2 ) is less sensitive 
to r\: variations of ij by a factor of three at the epoch 
t > 40 yr (see FigJ^J) only lead to a variation of the CR 
pressure P c /(po^, 2 ) ~ 1 — Os/c by a factor of 1.5. Since 
the CR pressure is dominated by the highest energy CRs 
one should expect a corresponding variation of the high 
energy gamma-ray emission. This is indeed seen from 
FigEJ for t > 40 yrs the expected gamma-ray flux at en- 
ergies e 7 — 1 — 10 TeV for low injection rate is lower by 
a factor ~ 1.5 compared with the case of high injection. 

Therefore, taking into account the uncertainties of all 
relevant parameter values we expect that the prediction 
for the TeV-emission is uncertain at best by a the factor 
of two (Fig 

5. SUMMARY 

A kinetic nonlinear model for CR acceleration in SNRs 
has been applied in detail to SN 1987A, in order to com- 
pare its results with observed properties. It is found that 
quite reasonable consistency with most of the observa- 
tional data can be achieved. 

The evidence for efficient CR production, leading to 
a strong modification of the shock, comes from the ra- 
dio synchrotron data. A proton injection rate of size 
rj « 10~ 3 is required to produce a significant shock mod- 
ification that leads to the steep spectrum of energetic 
electrons which fits the observed synchrotron spectrum 
very well. The condition is an extremely high down- 
stream magnetic field strength ~ 10 mG. Such a 
high field implies significant synchrotron losses of CR 
electrons emitting nonthermal X-rays. This makes the 
high frequency part of the synchrotron spectrum much 
softer and is consistent with the high-energy part of the 
X-ray spectrum, which is presumably of nonthermal ori- 
gin. Therefore the fit of the synchrotron spectrum gives a 
good estimate for the CR injection parameters and makes 
it possible to calculate the expected 7-ray flux in a spher- 
ically symmetric model. However, the basic tendency of 
the magnetic field vectors to form an Archimedean spiral 
configuration in the circumstellar wind interaction region 
requires a renormalization of the overall 7-ray flux. This 
a posteriori reduction by a factor / ro = 0.2 is chosen 
following the previous analyses of other SNRs. This in- 
troduces an acknowledged uncertainty in the predicted 
overall 7-ray flux. 

Since the SN shock interacts with the cool shell which is 
the densest part of the CSM, the 7r°-decay 7-ray spectral 
energy flux density at the current epoch is already quite 
high e 1 F 1 ss 4 x 10~ 13 erg cm~ 2 s _1 at energies e 7 = 
0.1 - 10 TeV. Depending upon the details of the CSM 
distribution during the next 15-20 yr the 7-ray flux is 
expected to increase roughly by a factor of two. 

The further temporal evolution of the gamma-ray emis- 
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sion corresponds to a secular decrease because the gas 
density in the unperturbed Red Supergiant Wind region 
decreases with radius oc r~ 2 . Since the particle injec- 
tion into the shock acceleration process can not be well 
determined for this later phase, the precise form of the 
emission decrease is not well known. 

The detection of 7-ray emission from SN 1987A would 
be a very important element in a consistent picture for 
this SNR. In particular, it would give evidence for ef- 
ficient CR production followed by strong magnetic field 
amplification for a core collapse supernova at a very early 
stage of evolution of its remnant. 



Accro, F., et al. 2010, A&A, 516, A62 
Aharonian, F., et al. 2006, A&A, 449, 223 
Bell, A. R. 2004, MNRAS, 353, 550 

Berezhko, E. G., Elshin, V. K., & Ksenofontov, L. T. 1996, JETP, 
82, 1 

Berezhko, E.G. & Volk, H.J. 2000, A&A, 357, 183 
Berezhko, E. G., & Ksenofontov, L. T. 2000, Astronomy Letters, 
26, 639 

Berezhko, E. G., Puhlhofer, G., Volk, H. J. 2003, A&A, 400, 971 
Berezhko, E. G., & Ksenofontov, L. T. 2006, ApJ,650, L59 (BK). 
Berezhko, E.G. 2005, Adv. Space Res., 35, 1031 
Berezhko, E.G. 2008, Adv. Space Res., 41, 429 
Berezhko, E.G., Ksenofontov, L.T., & Volk, H.J., 2009, A&A , 505, 
169 

Berezhko, E.G. & Volk, H.J. 2010, A&A, 511, A34 

Blondin, J.M. & Lundqvist, P. 1993, ApJ, 405, 337 

Chevalier, R. A. & Dwarkadas, V. V. 1995, ApJ, 452, L45 

Chevalier, R. A., & Fransson, C. 1987, Nature, 328, 44 

Crots, A.P.S., Kunkel, W.E. & Heathcote, S.R. 1995, ApJ, 438, 724 

Dorfi, E. A. 1991, A&A, 251, 597 

Drury, L.O., et al. 1999, A&A, 347, 370 

Ellison, D.C., et al. 2010, ApJ, 712, 287 

Enomoto, R., et al. 2007, ApJ, 671, 1939 

Gaensler, B. M., Manchester, R. N., Staveley-Smith, L., et al. 1997, 
ApJ, 479, 845 

Gorenstein, P., Hughes, J. P., Tucker, W. H. 1994, ApJ, 420, L25 
Hasinger, G., Ashcnbach, B. & Triimper, J. 1996, A&A, 312, L9 
Komin, N., et al. (H.E.S.S. Collaboration) 2010, at COSPAR 

Scientific Assembly, Bremen, E19-0096-10 (Poster, Nr. Thu-299) 
Ksenofontov, L. T., Volk, H. J., & Berezhko, E. G. 2010, ApJ, 714, 

1187 



This work has been supported in part by the Rus- 
sian Foundation for Basic Research (grants 09-02-12028, 
10-02-00154), the Council of the President of the Rus- 
sian Federation for Support of Leading Scientific Schools 
(project no. NSh-3526.2010.2), the Ministry of Edu- 
cation and Science (contract 02.740.11.0248). The au- 
thors thank Dr. V.N. Zirakashvili for pointing out a nu- 
merical inaccuracy in the calculation of the X-ray syn- 
chrotron spectra in an early version of this paper. EGB 
acknowledges the hospitality of the Max-Planck-Institut 
fur Kernphysik, where part of this work was carried out. 



Lucck, S. G. & Bell, A. R. 2000, MNRAS, 314, 65 
McCray, R. 1993, ARA&A, 31 175 

Murphy, T., Gaensler, B. M.& Chatterjee, S.2008, MNRAS,389, 
L23 

Ng, C.-Y., Gaensler, B. M., Stavclcy-Smith, L. et al. 2008, ApJ, 
684, 481 

Park, S., Burrows, D. N., Garmire, G.P. ct al. 2007, in AIP Conf. 

Proc. 937, Supernova 1987A: 20 Years After ed. S. Immler, K.W. 

Weiler, & R. McCray (Melville, NY: AIP), 43 
Potter, T. M., Staveley-Smith, L., Ng, C.-Y. et al. 2009, ApJ, 705, 

261 

Ptuskin, V. S. & Zirakashvili, V. N. 2003, A&A, 403, 1 
Ptuskin, V. S. & Zirakashvili, V. N. 2005, A&A, 429, 755 
Staveley-Smith, L., Manchester, R. N., Kesteven, M. J. ct al. 1992, 

Nature, 355, 147 
Turtle, A. J., Campbell- Wilson, D., Bunton, J. D., ct al. 1987, 

Nature, 327, 38 

Volk, H.J. 2004, in Frontiers of Cosmic Ray Science, Proc. 28th 
ICRC, Tsukuba, 8, 29 ff, ed. T, Kajita, Y. Asaoka, A. Kawachi, 
Y. Matsubara, & M. Sasaki (Tokyo, Japan, Universal Academy 
Press Inc.) 

Volk, H. J., Berezhko, E. G., & Ksenofontov, L. T. 2003, A&A409, 
563 

Volk, H. J., Berezhko, E. G., & Ksenofontov, L. T. 2005, A&A, 
433, 229 

Zanardo, G. et al. 2010, ApJ, 710, 1515 

Zhekov, S. A., Park, S., McCray, R., et al. 2010, MNRAS,407, 1157 
Zirakashvili, V. N. & Aharonian, F.A. 2010, ApJ, 708, 965 



